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Qs, . Abstract. 

We study, via numerical experiments, the role of bound states in the evolution of 
cosmic superstring networks, being composed by p F-strings, q D-strings and (p, q) 
bound states. We find robust evidence for scaling of all three components of the 
network, independently of initial conditions. The novelty of our numerical approach 
^-i, consists of having control over the initial abundance of bound states. This indeed 

<-j | allows us to identify the effect of bound states on the evolution of the network. Our 

studies also clearly show the existence of an additional energy loss mechanism, resulting 
to a lower overall string network energy, and thus scaling of the network. This new 
^ ' mechanism consists of the formation of bound states with an increasing length. 
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1. Introduction 

Cosmic superstrings can be formed as a result of brane annihilations in the context of 
brane- world scenario [U El E] • In the brane- world models, brane inflation takes place pi] 
while two branes move towards each other, and their annihilation releases the brane 
tension energy that heats up the universe to start the hot big bang era. Typically, 
strings of all sizes and types may be produced during the collision. Considering a IIB 
string theory in a (9+l)-dimensional space-time, interactions of Dirichlet (D) branes 
leads to the unwinding and subsequent evaporation of higher dimensionality branes, 
with the survival of three-dimensional (D3) branes embedded in a (9+l)-dimensional 
bulk — one of which could play the role of our universe — and Dl-branes (D-strings) [5]. 
Large Fundamental (F) strings and D-strings that survive the cosmological evolution 
become cosmic superstrings. They are of cosmological size and could play the role of 
cosmic strings [6J UJ, false vacuum remnants formed generically at the end of hybrid 
inflation within Grand Unified Theories [HI EJ - Cosmic superstrings have gained a lot of 
interest, particularly since it is believed that they may be observed in the sky, providing 
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both, a means of testing string theory and a hint for a physically motivated inflationary 
model (for a recent review, see e.g. Ref. |10j). 

Brane collisions lead also to the formation of bound states, (p, g)-strings, which 
are composites of p F-strings and q D-strings [HI [12]. The presence of stable bound 
states implies the existence of junctions, where two different types of string meet at 
a point and form a bound state leading away from that point. Thus, when cosmic 
superstrings of different types collide, they can not intercommute, instead they exchange 
partners and form a Injunction, as a consequence of charge conservation at the junction 
of colliding (p, g)-strings. The evolution of F, D strings and their bound states is a 
rather complicated problem, which necessitates both numerical as well as analytical 
investigations. Junctions may prevent the network from achieving a scaling solution, 
invalidating the cosmological model leading to their formation. In a number of studies, 
cosmic superstring evolution has been addressed via numerical experiments [13l [HI HHJ 
IT6l [T7[ [T9] . The formation of three-string junctions and kinematic constraints for their 
collisions have been also investigated analytically [201 |2T[ 122] . 

In what follows, we address the question of the effect of junctions in the evolution 
of a cosmic superstring network, being composed by three components: p F-strings, 
q D-strings and their (p, q) bound states. The results of our numerical investigations 
can be summarised as follows. Firstly, there is clear evidence for scaling of all three 
components of the network, independently of the chosen initial configurations. Secondly, 
the existence of bound states effects the evolution of the network. Thirdly, for (p, q) 
strings there is a supplementary energy loss mechanism, in addition to the chopping off 
of loops, and it is this new mechanism that allows the network to scale. More precisely, 
the additional energy loss mechanism is the formation of bound states, whose length 
increases, lowering the overall energy of the network. We note that in our simulations 
we can have control over the initial population of bound states, and this renders our 
novel results particularly important. 

2. The model 

We adopt the model of Ref. [17]. This is a simple field theory model of (p,q) bound 
states, in analogy to the Abelian Higgs model, which incorporates the main features of 
string theory. To represent the two different species of strings, the model includes two 
complex scalar fields, <fi and \i °f the Abelian Higgs model, coupled via a potential, 
so that bound states can be formed. The presence of stable bound states implies the 
existence of junctions, where two different types of string meet at a point and form a 
bound state leading away from that point. 

We distinguish two different cases of cosmic superstring networks. The nature 
of cosmic superstrings depends on the particular brane inflationary model leading to 
their formation. In the case that both species of cosmic strings are BPS, the model is 
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described by the action [T7| : 

S= J d 3 xdt[- \F 2 - \ (D,0) (D^Y - h (00* - 

-\ R2 -\ ( D »X) (D» X )* - (xx* - % 2 ) 2 ] , (1) 

where the covariant derivative is defined by 

D^cj) = d^cj) - ie x A^ , 

D^X = d^x ~ i^C^x ■ (2) 

For clarity in our discussion, we label the field as "Higgs" and the X field as 
"axion", even though both fields are Higgs-like. The scalars are coupled to the U(l) 
gauge fields and C^, with coupling constants e\ and e 2 and field strength tensors 
F\us — d fl A u — d u A^ and H^ v = d^C u — d u C^, respectively The scalar potentials are 
parametrised by the positive constants Xi,r]i and \2,V2, respectively. 

In the case that one species of string is non-BPS, we remove the second gauge 
field by setting e 2 = 0. In this way, this species of string is represented by the 
topological defect of a complex scalar field with a global U(l) symmetry. Such defects 
are characterised by the existence of long-range interactions [18J — as opposed to local 
strings in which all energy density is confined within the string, so that local strings 
have only gravitational interactions — implying different consequences for the evolution 
of the network. 

In flat space the classical equations of motion for the and X fields, obtained from 
the above action, Eq. ([1]), read [17] 

dpF 1 ™ =2eilm0*/r0, 

=2e 2 lm X *D u X , 
D^<j> = - 2 A x (0*0 - V l) - A 2 ( X *x -v!) 2 <P 

D,D» X = -2A 2 0*0( X *X-% 2 )X- (3) 

By calculating numerically the tension of the (p, q) strings, for a range of values of the 
charges, and then fitting the data to the square-root expression 



Hp,i) = /Wp 2 + ?7#s , ( 4 ) 

where /^f denotes the effective fundamental string tension after compactification and g s 
stands for the string coupling, we check [TTJ the validity of our model. We remind to 
the reader that Eq. (j4j) is exact only in the BPS limit. 

The potential is chosen such that bound states are energetically preferred over single 
strings. This can be seen from the form of the potential for the axion strings. At the 
core of the Higgs string the field vanishes, = 0, and therefore the potential of the 
axion vanishes as well. Therefore it is energetically favourable for the axion string to 
be located at the core of the Higgs string. In the true vacuum of the Higgs, we have 
|0| = 771, so the axion string behaves just like the Higgs string. 
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We employ this model to study, via numerical simulations, overall properties of the 
cosmic superstring network. In particular, we are interested in investigating the effect 
of junctions and their role in achieving scaling of the network. 

In the actual simulations we include the expansion of the universe by considering 
strings of fixed co- moving thickness, according to the model of Refs. [23j [2U [25]. In 
a Friedmann-Lemaitre- Robertson- Walker universe and considering the temporal gauge, 
A = 0, the equations of motion for the scalar fields take the form: 

d F 0i + 2 (1 - s) ^F 0l =2a 2s e 1 Im<J)*D u <J) + d J F 3i (5) 



d C 0i + 2(l-s)-C 0l =2a 2s e 2 luix f D v X + d j C 3i (6) 

2A 1 (0*0-r ?1 2 )+A 2 ( X *X-% 2 ) 2 ]0 (7) 
X " + 2-x' + A Ax = - 2A 2 a 2 *0*0 ( X *X ~ r,l) X (8) 



a 

a' 



b" + 2-0' + AA0 = - a 2s 
a 

y 

a 

where the time variable is the conformal time r: 



ds 2 = a 2 (r) (-dr 2 + dx 2 ) . (9) 

In our simulations we set s = 0, instead of the s = 1 required by the continuum-limit 
of the equations of motion (otherwise the defects will fall through the lattice); we also 
take, for simplicity, a' /a = 1/r. 

We calculate the correlation length using the expression: 

e = y|, (10) 

where V is the volume of the simulation box and L is the total length of the string 
network inside the box (following the same approach as for example in Ref. [IH]). We 
identify the location of the string by the value of the scalar fields. 



3. Results 



We present our main results for the scaling of (p, q) string networks, using simulations 
of the classical time evolution on lattices with three spatial dimensions. We consider 
two types of networks, namely local-global and local-local networks. The axion field has 
a global U(l) symmetry, and therefore the axion strings display long-range interactions. 

The lattice discretisation of the equations of motion is discussed in Ref. [17] . We 
choose all parameters of the Lagrangian to have natural values, 

771 = 1, Ai = 1, e x = 1, 772 = 1, A 2 = 1, e 2 = 1 . (11) 

The lattice spacing is 5x = 0.45, roughly half the characteristic length scale set by the 
masses of the scalar and gauge fields in the broken phase. The time-step is 5t = 10~ 2 . 
The simulations are carried out with boxes of volume 256 3 . 

For each of the two (local-global and local-local) networks we consider two types 
of initial conditions: one in which a large percentage of strings (~ 50 % of the string 
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length) is in bound states, and the other one with a small percentage of strings (~ 2 % 
of the string length) is in bound states. We are able to achieve the small bound state 
percentage configuration, by starting with a large percentage configuration and then 
shifting one of the fields by roughly half the initial correlation length. 

The type of behaviour we observe is that in most cases the string network does have 
a scaling behaviour, but the slope, 7, of the correlation length, £, versus time, defined 
by: 

£(r)= 7 r, (12) 

undergoes a discrete change as the network evolves. This change is most clearly visible 
in the case of a local-global string network, for both large and small amounts of bound 
states. The most dramatic change is that for the correlation length of the bound states 
themselves, provided the amount of bound states is small. 

For local-local networks, the change in correlation length is not detectable if there 
is a large amount of bound states, while if there is a small amount of bound states we 
find a case where the bound states themselves do not exhibit scaling behaviour. This 
violation of scaling by the bound states can be explained by taking into account that the 
formation of bound states gives an additional possibility for a network of (p, q) strings to 
decrease its energy. Therefore, the length of the bound states increases as a function of 
time, violating the scaling behaviour which requires the total length to decrease. Once 
the length of the bound states have reach the maximum value, scaling resumes. This 
violation of scaling by the bound states does not lead to violation of scaling for the 
entire network, as it can be observed from the evolution of the Higgs and axion strings. 

In Fig. [H we show the string correlation length for the Higgs and axion fields, as well 
as for their bound states, as a function of time. The initial configuration is a local-global 
network with a large amount of bound states. Clearly, there is convincing evidence for 
scaling of the three components of the network. This scaling is characterised with a 
distinct change of the correlation length slope during the network evolution. 




Figure 1. The Higgs (left), axion (middle), and bound state (right) string correlation 
length as a function of time. The network is a local-global one with a large amount of 
bound states. The data and linear fits for the two regimes are shown. 



In Fig. [2] we show the string correlation length for the Higgs and axion fields, as 
well as for their bound states, as a function of time. The initial configuration is a 



Dynamics of F/D networks: the role of bound states 



6 



local-global network with a small amount of bound states. Again, scaling behaviour is 
apparent for all three components of the local-global network. The difference between 
the case plotted in Fig. [2] and the one plotted in Fig. CEJ is the percentage of initial 
bound states. As one can easily observe from Fig. [2j if the initial state of a local-global 
network has a small amount of bound states, then the change of the correlation length 
slope during the evolution of the network is more acute. 




ao Time i?o 140 60 no Time i?o ho w bo Time 120 -so 



Figure 2. The Higgs (left), axion (middle), and bound state (right) string correlation 
length as a function of time. The network is a local-global one with a small amount of 
bound states. The data and linear fits for the two regimes are shown. 



We then draw the corresponding plots for local-local networks, with a large and a 
small amount of bound states in Fig. [3] and Fig. H] , respectively. We see, from Fig. [31 
that for a local-local network with a large amount of bound states, all three components 
scale (as for the local-global case), whereas the change in correlation lengths during 
evolution is not really apparent (as opposed with the local-global case). 
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Figure 3. The Higgs (left), axion (middle), and bound state (right) string correlation 
length as a function of time. The network is a local-local one with a large amount of 
bound states. The data and linear fits for the two regimes are shown. 



In Fig. HI we plot the correlation length of the Higgs, the axion and their bound 
states for a local-local network with a small amount of bound states. Here we explicitly 
see that there is a violation of the scaling behaviour for the bound states, as a result 
of the increase of their length with time. As we have already mentioned earlier, this 
violation does not affect the entire network, which indeed scales. 

In Fig. [5] we show, for comparison, the correlation length of a single species of 
string, using the same initial data as the one used for the (first) Higgs of the (p, q) 
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Figure 4. The Higgs (left), axion (middle), and bound state (right) string correlation 
length as a function of time. The network is a local-local one with a small amount of 
bound states. The data and linear fits for the two regimes are shown. 



strings. Clearly, there is no change in the slope of the correlation length as a function 
of time. This is an important check to assure that finite size effects do not influence the 
results of the simulations. 



60 so Time 120 140 



Figure 5. The correlation length for a single species of cosmic string. There is no 
change in the slope of the correlation length as a function of time. 



4. Conclusions 

In the context of brane-world models, brane collisions lead to the formation of 
large Fundamental (F) strings and/or Dl-branes (D-strings). Those who survive the 
cosmological evolution become cosmic superstrings, playing the role of their solitonic 
analogues. By observing strings in the sky, we may be able to test for the first (and 
maybe only) time string theory. To know the observational consequences of cosmic 
superstrings, it is essential to master the properties and evolution of such networks. 

Performing numerical experiments, we have studied the dynamics and overall 
properties of superstring networks. More precisely, we have performed field theory 
simulations of p F- and q D-strings, and their (p, q) bound states. We have investigated 
the effect of bound states and in particular the approach to scaling, which is crucial for 
the cosmological consequences of cosmic superstring networks. Defects which do not 
scale are cosmologically undesired, since they may over-close the universe, leading to 
the old monopole problem. 
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Our studies have clearly shown that the three components of the network scale, 
independently of the chosen initial conditions. In addition, by having control over the 
initial abundance of bound states, we were able to identify the effect of bound states 
on the overall network evolution. Finally, we have found that there is an additional 
energy loss mechanism, beyond the chopping off loops. This new mechanism consists on 
the formation of bound states, resulting to a lower overall energy of the network, thus 
leading to scaling. 
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